Neutron scattering from a coordination polymer quantum paramagnet 
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Inelastic neutron scattering measurements are reported for a powder sample of the spin- 1/2 quan- 
tum paramagnet Cu(QuinoxaIine)Br2. Magnetic neutron scattering is identified above an energy 
gap of 1.9 meV. Analysis of the sharp spectral maximum at the onset indicates that the material is 
magnetically quasi-one-dimensional. Consideration of the wave vector dependence of the scattering 
and polymeric structure further identifies the material as a two-legged spin- 1/2 ladder. Detailed 
comparison of the data to various models of magnetism in this material based on the single mode 
approximation and the continuous unitary transformation are presented. The latter theory provides 
an excellent account of the data with leg exchange Jy = 2.0 meV and rung exchange J± — 3.3 meV. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 



I. INTRODUCTION 

Progress towards understanding the cooperative quan- 
tum physics of one dimensional systems is frequently 
gated by the availability of suitable model systems for 
experiments. While oxides offer the possibility of car- 
rier doping, coordination polymer magnets have energy 
scales that are well suited for high magnetic field experi- 
ments that traverse phase boundaries. Such experiments 
have been important for elucidating the fermionic na- 
ture of low energy excitations in the uniform spin-1/2 
chain^ and their remarkable transition to a gapped soli- 
tonic spectrum upon application of an effective staggered 
field.^ 

A natural next step in the experimental exploration of 
one dimensional spin-1/2 systems is to examine magnetic 
excitations as a function of field and temperature in two 
legged spin-1/2 ladders. The Hamiltonian of this system 
has the form 



H = Y.JA,^■S2,i (1) 

i 

+ • Si_,;+i + S2^i • S2^i+l), 

where i indexes rungs and 1, 2 each leg as illustrated 
in Fig. n] The Lieb-Schultz-Mattis theorem^ allows 
for a finite gap in the excitation spectrum and indeed 
analytical;^i^iSiL2i2iiS numerical^ and experimental^^ 
work has confirmed that a spin gap does exist. A key 
experiment is to close this spin gap with an external 
field and explore the predicted quantum critical high field 
phaseiJ^ A rich zero field excitation spectrum that should 
feature both bound states and continua also remains to 
be examined experimentally>iS*H 



Unfortunately as of now no experimental model system 
has been identified to enable neutron scattering experi- 
ment in the high field phase. (La, Sr, Ca)i4Cu2404i— 
and SrCu203ii do appear to be spin ladders but the 
high energy scales for these systems render them unsuit- 
able for the proposed experiments. The energy scales for 
Cu2(l, 4 - diazacycloheptane)2Cl4 (CuHpCl)iS-iSiiiii2iia 
and (VO)2P20 , are appropriate and they were previ- 
ously identified as spin-1/2 ladders on the basis of specific 
heat and magnetization measurements. However, further 
investigations with neutron scattering revealed that the 
exchange interactions in these materials do not corre- 
spond to spin ladders despite spin ladder structural mo- 
tifs. In the case of CuHpCl there were early indications 
from magnetization measurements that the material is 
not one dimensional and subsequent neutron scattering 
experiments showed that the group velocity for magnetic 
excitation has no substantial anisotropy.'^^ In addition, 
the molecular bond, which was thought to be the rung 
of the ladder, assumes a frustrated configuration in the 
ground state. Hence CuHpCl is now thought to have 
a complex frustration induced singlet ground state,— 
(VO)2P207 on the other hand is quasi-one-dimensional 
but rather than being a spin ladder, neutron scattering 
experiments showed that it is an alternating spin chain 
with the chain direction perpendicular to the putative 
ladder directionj^^i^^i^^ 

These false starts illustrate that bulk measurements, 
while sensitive to the magnetic density of states, can 
not provide definite conclusions regarding the geometry 
of interacting spin system with a singlet ground state. 
With the goal of identifying spin ladder systems for high 
field experiments, we have initiated a project to synthe- 
size potential coordination polymer spin ladders and sub- 
ject them to inelastic neutron scattering experiments for 
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FIG. 1: Schematic of a two-leg spin ladder system, indicai 
ing the rung interaction J± and the rail interaction Jy. Th 
indexing of spin sites used in Eq. ^ is also shown. 



verification. Cu(Quinoxaline)Br2 is a first out of se\ 
eral such candidate materials!^ Susceptibility measure 
ments on this materialSi are well described by a spir 
ladder model with antiferromagnetic exchange constant 
J|| = 2.61 meV and J± — 3.02 meV. Magnetizatio 
measurementsSi show a phase transition from a nor 
magnetic singlet ground state to a magnetized state nea 
a critical field He = 14 T, indicating the presence c 
an energy gap A = giisHc ~ 1.7 meV. Furthermor 
the square root singular onset of magnetization abo'v 
the critical field suggests that the system is magneticall 
one-dimensional. 

In the present paper, we report inelastic neutron seal 
tering experiment on a deuterated powder sample c 
Cu(Quinoxaline)Br2 to determine whether this materif 
is indeed a spin ladder and to establish the strength of th 
relevant exchange interactions. The comparison of wav 
vector and energy dependent magnetic neutron scatter- 
ing data to comprehensive model calculations carried out 
as a function of rung and leg exchange constants strongly 
indicate that Cu(Quinoxaline)Br2 is a coordination poly- 
mer spin ladder, with weak inter-ladder interactions. Fi- 
nal verification will require inelastic scattering experi- 
ments from single crystals. 

The crystal structure of Cu(Quinoxaline)Br2 
(Cu(C8H6N2)Br2) is monoclinic with space group C2/m 
and lattice constants a = 13.1745(15) A, b = 6.9293(8) 
A, c = 10.3564(12) A, and /3 = 107.699(2)°^^ Quinox- 
aline (C8H6N2) has a tendency towards polymeric 
structures where quinoxaline forms a bridge between 
metal atomsiSS In Cu(Quinoxaline)Br2, Cu2Br4 dimers 
are linked to adjacent dimers by bridging quinoxaline 
molecules along the monoclinic b axis as shown in 
Fig. OTa). The inter-dimer spacing along the b-axis is 
6.929 A. Hydrogen bond mediated interactions may ex- 
ist between molecular units displaced by di = (a±b)/2 
and d2 = c (See Fig.|2Jl. If there were sufficiently strong, 
the material could be a two-dimensional bi-layer system 
as in BaCuSi206r^ an alternating spin chain as in 
Cu(N03)2 • 2.5D20r^ or a three dimensional system of 
interacting spin pairs as in TlCuCla^ and KCuCl3>2Si^ 




FIG. 2: Structure of Cu(Quinoxaline)Br2 projected along 
four different directions to illustrate potential interact- 
ing spin models. (a) The molecular ladder structure of 
Cu(Quinoxaline)Br^-, which the ladder extends along the 
monoclinic b axis. Interactions along the (a ± b)/2 and c 
directions could lead to two dimensional dimer models in the 
a-b or b-c planes (frame (b) and (c)) or a three dimensional 
dimer model (frame (d)). Color coding is as follows: Pink:Cu, 
GreeniBr, Blue:C, Red: H, and Grey:N. 



II. EXPERIMENTAL TECHNIQUES 

A. Sample and Neutron Instrumentation 

The neutron scattering sample consisted of 27 grams of 
deuterated powder in an annular aluminum can with in- 
ner diameter 2.42 cm, outer diameter 2.92 cm, and height 
10 cm. The powder was obtained through precipitation 
of a 1:1 molar solution of anhydrous CuBr2 and deuter- 
ated d-4 quinoxaline dissolved in a small amount of 95% 
ethanol. Prompt gamma neutron activation analysis and 
NMR measurements showed that 58% of the hydrogen 
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sites in our sample were occupied by deuterium. This is 
less than the 2 /3 deuteration level of the d-4 quinoxaline 
starting material indicating some hydrogen/deuterium 
exchange with the solvent. 

Inelastic neutron scattering measurements were per- 
formed using the disk chopper time-of-flight spectrom- 
eter (DCS)^^ at the National Institute of Standards 
and Technology (NIST) Center for Neutron Research in 
Gaithersburg, Maryland. A disk chopper system was 
used to select a 167 Hz pulsed neutron beam with an 
energy of 4.87 meV and a pulse width of 79 /iS from the 
NCNR cold neutron source. The 0.65 steradian detection 
system probed inelastic scattering with energy transfer 
—4.17 meV < hio < 3.84 meV and momentum transfer 
0.13 A-i < Q < 3.40 A-i. The full width at half maxi- 
mum (FWHM) elastic energy resolution was Shui ~ 0.18 
meV. 



B. Neutron Scattering Data Analysis 

The remnant hydrogen in the sample gives rise to a 
scattering cross section of 201 barn per formula unit com- 
pared to a total magnetic scattering cross section from 
Cu^"^ of 0.146 barn per formula unit. There are several 
distinct challenges associated with the strong incoherent 
nuclear scattering. (1) The strong elastic scattering cross 
section reveals tails of the energy resolution function that 
can dwarf the magnetic signal well beyond the FWHM 
of the resolution function. (2) Phonon scattering and 
multiple scattering processes involving phonons produce 
a strong background even at low temperatures that must 
be determined accurately in order to isolate the mag- 
netic scattering. (3) Both magnetic and phonon scat- 
tering are frequently preceded or followed by incoherent 
elastic nuclear scattering. Because incoherent scattering 
is approximately wave vector independent this leads to 
a Q-averaged contribution of both types of scattering at 
any wave vector. We have developed accurate techniques 
for dealing with each of these issues and because they are 
essential for this research they are described in the fol- 
lowing sections. 



1. Temperature independent background 

The elastic background from the tails of the resolution 
function can be determined by utilizing the fact that in- 
elastic scattering obeys the principle of detailed balance 
whereas elastic scattering at sufficiently low tempera- 
tures, to a good approximation, can be approximated 
as being temperature independent. Under these circum- 
stances the raw measured count rate, Ir{uj,T), for any 
specific value of wave vector transfer satisfies: 

/,(c^,r) = BiLo)+Iiio,T). (2) 
IrHu^lT) = Bi-\u;\) + Ii\u;lT)eM-fiH/kT)i3) 



Here B{u;) is the temperature independent background 
and I(uj,T) represents all inelastic scattering processes 
that satisfy detailed balance at the temperature T. An 
underlying assumption here is that resolution effects can 
be neglected for the inelastic scattering. By measuring 
data at two different temperatures ri=1.4 K and r2=60 
K and using equations Q and Q it is possible to extract 
values for B{~\uj\), I{uj,Ti) and/(w,T2) from the 

corresponding four equations. 



2. Subtracting phonon scattering 

The above procedure isolates inelastic scattering from 
temperature independent elastic scattering and detector 
dark current. However, this scattering intensity still has 
contributions from magnetic as well as phonon scatter- 
ing. To remove the low temperature contribution from 
phonon scattering, Ip{u),Ti) , and hence isolate low tem- 
perature magnetic scattering, Im{i-^,Ti), we use the fact 
that phonon scattering in hydrogenous systems dwarfs 
magnetic scattering for ksT >> huj. Hence we make the 
following approximation for kBT2 > tko: 

/(W,T2) = /p(c^,T2) +/,„(C^,T2) « Ip{uJ,T2). (4) 

We further assume that one phonon scattering events 
dominate over multi-phonon events and we neglect the 
temperature dependence of the Debye Waller factor such 
thaliS 



/p(w,Ti) 



1 — exp(— 

1 — exp(— /5i?icj) 



(5) 



The final expression for the magnetic inelastic scatter- 
ing as deduced from the measurements at Ti = 1.5 K 
and T2 = 60 K is therefore 

T I rr,^ Tf ^ 1 ~ eXp(-^2fiw) w „ ^ 

1 — exp(— /Jiftw) 



3. Multiple Scattering 

Multiple scattering is difficult to avoid when probing 
weak magnetic scattering in the presence of strong in- 
coherent scattering. Here we describe how to account 
for the dominant double scattering process that involves 
incoherent elastic scattering preceded or followed by a 
weaker inelastic process. For simplicity we assume that 
the incoherent elastic scattering is isotropic and inde- 
pendent of neutron energy and we neglect anisotropy 
introduced by the sample geometry.^^ Under those cir- 
cumstances the incoherent scattering event in a double 
scattering process effectively randomizes the direction of 
scattering for the preceding or subsequent inelastic scat- 
tering: 



I'{Q,u;) = TIiQ,Lu) 



(7) 
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FIG. 3: Scaled wave vector dependence of inelastic neutron 
scattering from Cu(Quinoxaline)Br2 at T=60 K and for vari- 
ous values of energy transfer. The scaling behavior is consis- 
tent with double scattering by phonons and incoherent elastic 
nuclear scattering as described by Eq. ^ and ||HJ. 



Here ki and kf are the incident and scattered wave vec- 
tors respectively. The integral implements the average 
over scattering directions for the inelastic process. T 
is an "effective" sample transmission that can be deter- 
mined through analysis of phonon scattering. Approxi- 
mating the intensity of single event phonon scattering for 
a thin sample as Ip(Q,u!) = f{u!)Q^ we find that 



/;(Q,u;) = /H(fc2 



1,2 



+ (1-^))- (8) 



Fig. 13 shows phonon scattering intensity at T = 60 K 
as a function of Q^/{kf + k'j). Data for several values 
of energy transfer, huj, scale in agreement with Eq. © 
and yield an effective value of T = 0.40. For compari- 
son the average transmission through a spherical version 
of our sample as calculated considering only incoherent 
scattering is 0.30. In the following we shall use Eq. [7| 
with T ~ 0.40 when comparing the measured inelastic 
magnetic scattering, Im(uj,T), to theoretical models of 
spin dynamics in Cu(Quinoxaline)Br2. 



4- Normalized magnetic scattering 

Anisotropic self shielding effects associated with the 
annular sample geometry were taken into account using 
a numerical integration technique^ Absolute normaliza- 
tion of the data was subsequently achieved using elastic 
incoherent scattering from the sample, duly considering 
the hydrogen/deuterium ratio determined through neu- 
tron activation and NMR analysisiA^ The technique is 
estimated to be accurate to within 20%. The normal- 
ized magnetic scattering intensity, Im, is related to the 
resolution smeared dynamic spin correlation function as 



Im{Q,uj) = 2 J dQ'hduj"RQ^{Q-Q\uo-uo') 



x||F(Q')p5(Q',c.') 



(9) 



Here TZq^{Q — Q' , oj — to') is a unity normalized reso- 
lution function that is peaked on the scale of the FWHM 
resolution for Q k, Q' and fiio k, huj' , g ~ 2.12 is the 
Lande g-factor for Cu^+ in Cu(Quinoxaline)Br2.^^ ^'(Q) 
is the anisotropic magnetic form factor appropriate for 
a hole in the 3d^2_y2 orbitali^ Based on bond distances 
and their coplanarity it appears that the two short Cu-Br 
bonds and the two bonds to the neighboring quinoxaline 
molecules define the x-y plane for the 3dj.2_y2 copper 
orbitals. We used this assumption to model the effects 
of the anisotropic magnetic form factor. The spherical 
average of the dynamic spin correlation function is 



(10) 



and the dynamic spin correlation function proper is 

-E 



1 

27rfi! 



rr' 



< ^rW^^'IO) > e~*^ <^-^'\ (11) 



Unless specifically stated, we use the notation defined 
in S. M. Lovesey's book on thermal neutron scatteringi^ 



III. EXPERIMENTAL RESULTS 

Fig. ^ (a-b) shows raw Q-integrated data at T=60 K 
and T=1.4 K. The range covered is 0.5 A'^ < Q < 1.0 
A~^. The solid symbols in Fig.QJb) show B{iS) as deter- 
mined by Eq. Q and Eq. |(2Jl. It is clear from this anal- 
ysis that the strong incoherent elastic scattering from 
hydrogen produces intense tails of scattering intensity 
well beyond the 0.18 meV FWHM of the resolution func- 
tion. The increase of the count rate beyond 3 meV (see 
Fig.^Ja)) is attributable to the same effect through time 
of flight frame overlap. 

A small peak is visible in the raw low temperature 
data (Fig.^Jb)) at approximately 2 meV. Following back- 
ground subtraction, Fig.^d) clearly shows that this peak 
marks the onset of a continuum of scattering in the wave 
vector integrated spectrum. Fig. 2Jc) shows that this 
continuum is not present at T=60 K where as expected, 
thermally activated phonon scattering dominates. The 
solid line in Fig. EJd) shows the projected contribution 
from phonon scattering at 1.4 K as derived from the 
data in Fig. 0Jc) using Eq. Q. While incoherent in- 
elastic phonon scattering produces a background contri- 
bution that increases with hio^ it is apparent that the 
bounded continuum above a 2 meV spectral gap is 



5 




-3-2-10 1 2 3 
flcj (meV) 

FIG. 4: The energy dependence of Q-integrated neutron scat- 
tering from Cu(Quinoxaline)Br2 at (a) T=60 K and (b) T=1.4 
K before detailed balance correction. In both frames solid 
symbols show the temperature independent background de- 
termined as described in section U.B,1. Background sub- 
tracted data at (c) T=60 K and (d) T=1.4 K data. The 
solid line is phonon scattering contribution at T=1.4 K as 
determined by Eq. 0. 



associated with inelastic magnetic scattering. This spin- 
gap is consistent with the activation energy derived from 
bulk measurements^ 

Having verified that magnetic scattering can be iso- 
lated from incoherent nuclear scattering through the 
techniques described in section lTl Bl the analysis can now 
be replicated as a function of wave vector transfer. The 
results are summarized in Figs. [S] and El where the data 
have furthermore been normalized to report values for 
as defined by Eq. 

Fig. 13 shows the normalized scattering intensity 
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FIG. 5: Normalized powder inelastic neutron scattering inten- 
sity for Cu(Quinoxaline)Br2 at (a) T=60 K and (b) T=1.4 K 
after detailed balance correction as described in section U.B,1. 



I{Q,uj) for Cu(Quinoxaline)Br2 at T=1.4 K and T=60 
K after subtracting the temperature independent back- 
ground determined through Eq. |2Jl and Eq. At 
T=60 K the scattering increases with wave vector trans- 
fer as expected for phonon scattering. Further analy- 
sis of the Q— dependence of these data was provided in 
section IIIBI 3. At low temperatures there is a clear on- 
set of scattering for hui > 2 meV and this is consistent 
with magnetic neutron scattering from a quantum para- 
magnet. However, a growth in the scattering intensity 
towards higher Q confirms the conclusions from the pre- 
vious spectral analysis that even the low T scattering 
contains significant contributions from phonon scatter- 
ing. 

Fig.|Hfa) shows the phonon contribution to low temper- 
ature scattering derived from the T=60 K data based on 
Eq. |SJ|. Subtracting these data from the data in Fig.lSJb) 
concludes the process of isolating the Q and huj depen- 
dent magnetic scattering. From these data, which are 
displayed in Fig.|^b), it is apparent that there is a sharp 
onset of magnetic scattering for hu; ~ 1.9 meV with a 
finite-Q maximum. In the following we shall show that 
these are features of a quasi-one-dimensional antiferro- 
magnet with a singlet ground state and that the data 
furthermore are consistent with a spin ladder extending 
along the b-axis. 
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FIG. 6; Normalized inelastic scattering intensity for 
Cu(Quinoxaline)Br2 (a) Phonon scattering intensity at T—1.4 
K. (b) Magnetic scattering intensity at T=1.4 K after sub- 
tracting the phonon contribution as described in section 

n.B,i. 



IV. ANALYSIS AND DISCUSSION 

A. Single Mode Approximation 

For most but not all^ quantum magnets, the domi- 
nant spectral weight resides in a resonant mode. When 
this excitation is between a singlet ground state and an 
excited triplet without magnetic long range order, it is 
called a triplon.**^ We prefer to reserve the common term 
magnon for systems with magnetic long range order. The 
"single mode approximation" (SMA) provides an excel- 
lent account of the dynamic spin correlation function for 
a dominant triplon. The starting point is the following 
low temperature approximation for iS(Q,ti;), which con- 
centrates all spectral weight in a resonant mode. 

S"^{Q, u;) = S{Q)6{nuj - e(Q))5,0. (12) 

The first moment sum-rule^ links the equal time cor- 
relation function, iS(Q), to the dispersion relation as fol- 
lows 

S{Q) = -^^[Jd(Si,,-S2,,)(l-cos(Q-d))](13) 

Here we make the simplifying approximation that the 
molecular bond with displacement vector d dominates 



the structure factor. Within this framework, we exam- 
ine different models of magnetism in Cu(Quinoxaline)Br2 
by varying parameters in the following phenomenological 
dispersion relation that allows for inter-molecular inter- 
actions along the (a ± b)/2 and c directions: 

e(Q) = eb(fc) + 2Shfc cos 27rft, ■ cos 27r/c (14) 
+ Bi cos{2Trl) (15) 

Here Q = ha* + kh* + Ic* . For models that approach 
the spin ladder limit we use the following strong rung 
coupling perturbation expansion,*^ 

e-h{k) = J±[l + — cos27rfc+i(i)2(3-cos47rfc) 



8 Vi 



J± 4Vi- 
(2 cos 27rfc + 2 cos 47rfc — cos 6nk — 3) 



-o( 



J|| 



(16) 



For models where the dispersion along b is not dominant 
we use the simpler form 

eb(fc) So + SfcCos27rfc (17) 

The spherical average of 4S(Q,a;) (Eq. (|10|l ') was in 
general calculated numerically except for truly one- 
dimensional models where we used the analytical 
resullii^. 



B. One-Triplon CUT Model 

For an ab-initio description of the dynamic spin 
correlation function for the two-leg ladder, we use a 
particle conserving Continuous Unitary Transformation 
(CUT).^?!^?!^^!^^ The CUT is realized in a perturbative 
fashion about the limit of isolated rung dimers and can 
yield the spectral weight and dispersion of the triplon 
elementary excitations, as well as the multi-triplon con- 
tributions. Here we are interested in the low-energy part 
of the dynamical structure factor for x = J\\/ J_l < 1- In 
this regime the dynamical structure factor is dominated 
by the one-triplon part. Two-triplon contributions are 
located at higher energies and carry only a small part 
of the total spectral weight (10% for x = 0.6). It is 
thus justified to restrict the discussion to the one-triplon 
part of the dynamical structure factor, which is charac- 
terized by the one-triplon dispersion ecuT(Q) and the 
one-triplon spectral weight a^(Q). The high-order series 
for both quantities are expressed in terms of an internal 
parameter of the system, the energy gap A, instead of 
X, and quantitative results for x < 1 are obtained by 
extrapolationi^Si^ 

The dynamic spin correlation function at T = and 
for w > is given by 

ocut('7,w) = — 



huj - ecuT(g) + 
TT {hhj ~ ecuT(g))^ + 5"^' 



1 



(18) 
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Here 5 = 0.02 meV was used to avoid divergences with- 
out introducing broadening beyond the experimental res- 
olution. The spherical average of the dynamic spin cor- 
relation function is obtained as follows: 



5, 



1 



CUT 



sin 9d0 



27r 



(19) 



X 25cut(Q||6,^)(1 - cos(d(3^ sin(/))). 

The integral over can be computed analytically lead- 
ing to 



(20) 



X [l-Jo(-^yW^)], 

where y — Qb cos 9 and Jo{x) is the zeroth order Bessel 
function. From the structure of the material, we have 
d = 3.75 1 and b = 6.93 1. 



C. Comparison between theory and data 

We now compare each of these models of magnetism 
in Cu(Quinoxaline)Br2 to the observed magnetic scatter- 
ing intensity. Multiple scattering involving inelastic mag- 
netic scattering followed or preceded by incoherent elas- 
tic nuclear scattering was added to the calculated single 
event scattering as described by Eq. (TJ with the experi- 
mentally determined effective transmission T = 0.4. We 
also allowed for an overall additive constant to account 
for any discrepancies in the background subtraction. 



1. Q and E cuts 

For an overview, Fig.|Z|shows the energy dependence of 
the magnetic scattering averaged over wave- vectors from 
0.5 A-i to 1.0 A-i as follows: 



f Q^IiQ,uj)dQ 
J Q^dQ ' 



(21) 



These data are a measure of the magnetic density of 
states and are sensitive to the dimensionality of the spin 
system. The dashed line in Fig. [7| shows a fit to the 
one dimensional CUT spectrum. The peak in the data 
is broader than this truly ID model. The dotted line 
shows the spectrum for a two dimensional model where 
Bq = 3.82 meV and Bhk = 1-02 meV, which produce 
a gap in the spectrum that is consistent with the data. 
This 2D model does not produce a peak at the spec- 
tral onset as observed in the experiment. A quasi-one- 
dimensional model (solid line in Fig. [7^ with Bq = 3.82 
meV, Bk = 1.84 meV, Bhk = 0.08 meV and Bi = 0.14 
meV on the other hand yields an acceptable fit. From 
this we may conclude that Cu(Quinoxaline)Br2 is mag- 
netically quasi one-dimensional. 
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FIG. 7: Energy dependence of the normalized magnetic scat- 
tering intensity for Cu(quinoxaline)Br2 averaged over wave 
vectors 0.5 < Q <1.0 A~^. Open circles are experimen- 
tal data. The dashed line indicates the one-triplon dynamic 
spin correlation function from the Continuous Unitary Trans- 
formation theory. The dotted line shows the spectrum for a 
two dimensional model with Bo ~ 3.82 meV and Bhk = 1.02 
meV. The solid line shows the spectrum for a quasi-one- 
dimensional ladder model with Bq = 3.82 meV, B^ = 1.84 
meV, Bhk = 0.08 meV, and Bi = 0.14 meV. AU models were 
convolved with the experimental resolution function and an 
overall constant was fit to account for any discrepancies in 
the background subtraction procedure. 
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FIG. 8: Wave-vector dependence of normalized magnetic scat- 
tering intensity for Cu(quinoxaline)Br2 at 1.4 K, integrated 
over energy transfers 1.85 meV < hw < 2.65 meV. Open cir- 
cles are experimental data. The solid and dash lines indicate 
the dynamic spin correlation function of the one-triplon CUT 
and SMA ID alternating spin chain models, respectively, con- 
volved with the experimental resolution function. Multiple 
scattering was added to the models as calculated from Eq. ((TJ 
using the experimentally determined value of T = 0.40. Note 
that data for Q < 0.5 does not include the full energy 
range due to the kinematic limitations of the experiment (see 
Fig.inj. 
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To establish the relevant one dimensional model, we 
examine the wave- vector dependence of magnetic scat- 
tering averaged over the energy range from 1.85 meV to 
2.65 meV (Fig.|SJ|. I{Q) is sensitive to spatial spin corre- 
lations and in particular the intra-dimer spin spacing, d. 
The data show a rounded maximum at Qo — 0.8 A"^, in- 
dicating that the spins forming the dominant correlated 
spin pairs are separated by d w tt/Qq = 4 A. For com- 
parison the nearest neighbor Cu-Cu separation within the 
Cu2Br4 molecular unit is 3.75 A. The solid line in Fig.|Sl 
shows the CUT ladder model prediction. The excellent 
agreement with the data is strong support for the ladder 
model. Conversely the c-axis alternating chain model for 
one-dimensional magnetism in Cu(Quinoxaline)Br2 pro- 
duces the dashed line which is clearly inconsistent with 
the data. 



2. Intensity contour maps and global fits 

We performed global fits of the dynamic spin correla- 
tion function associated with the models described in sec- 
tion IV to the Q and uj dependent magnetic neutron scat- 
tering from Cu(Quinoxaline)Br2. The results are sum- 
marized in Figs. 13- A least-square fit of the ladder 
single-mode approximation gave (Si^; • §2,^) — —0.37(1), 
(Si^i-Si.i+i) = —0.04(1). The former value is close to the 
expectation value of -3/8 for a spin dimer singlet ground 
state. We also found J\\ = 1.98(4) meV, J_l = 3.05(1) 
meV. The ratio J\\/J± ~ 0.65 indicates that the pertur- 
bative expression for the dispersion relation in Eq. H16(l 
is at the limit of applicability. Allowing for inter-ladder 
coupling through finite values for Bhk and Bi reduced 
the from 1.3 to 1.1 for Bhk = 0.08 meV and Bi = 0.14 
meV, which indicates that there are weak interactions 
between ladders. Fig. lllf b') shows the best fit which cor- 
responds to 5o = 3.82 meV, Bu = 1.84 meV, Bhk = 0.08 
meV, Bi = 0.14 meV, and Jd(Si,, • S2,,;) = -1.0 mcV. 

As opposed to the strong coupling perturbation ex- 
pansion, the CUT can provide results for S{CI,ljj) for 
J_L J|| . Fits were undertaken for fixed exchange ra- 
tios X = J\\/Ji_ ranging from 0.5 to 1. For each value 
of X, the only adjustable parameters was an overall scale 
factor and the constant background. Jj_ was selected so 
as to produce the observed value of A for each value of x 
as listed in Tabled Fig.^shows the x dependent smallest 
possible value for the goodness of fit criterium. The 
best fit was obtained for x = 0.6. The corresponding 
values of the exchanges constants are J|| = 1.98 meV 
and J_L = 3.3 meV. The entire powder averaged dy- 
namic spin correlation function for this model is shown 
in Fig. llOf a). Fig. llir a) compares the CUT model calcu- 
lations for X — 0.6 to the experimental data. Fig. Iior b) 
shows that the strong coupling dispersion relation ob- 
tained from the SMA fits is in fact also in good agreement 
with that derived from the CUT. 

Whether the corresponding scattering cross section is 
approximated by the CUT or the SMA, the ladder model 



1.4 




FIG. 9: = — 5^ {1°^^ — Ij^)^ /o-^ versus the ratio of rail 

i 

to rung exchange in the comparison of magnetic neutron scat- 
tering data for Cu(quinoxaline)Br2 to the Continuous Unitary 
Transformation theory of one-triplon excitations in a b-axis 
spin ladder. 



TABLE I: Values of ladder exchange constants for which the 
one triplon contribution to S{Q,u!) was calculated for com- 
parison to the data 



J|l (meV) 


J± (meV) 


X = J\\/J± 


x' 


1.56 


3.12 


0.5 


1.39 


1.77 


3.21 


0.55 


1.34 


1.98 


3.3 


0.6 


1.30 


2.22 


3.41 


0.65 


1.37 


2.46 


3.52 


0.7 


1.45 


2.96 


3.70 


0.8 


1.49 


3.29 


3.80 


0.866 


1.48 


3.47 


3.85 


0.9 


1.48 


3.98 


3.98 


1.0 


1.46 



clearly provides an excellent account of the data. The 
virtue of CUT in this context is that apart from the ex- 
change constants and an overall scale factor that in prin- 
ciple can be determined through the total moment sum 
rule, there are no adjustable parameters. 



V. CONCLUSION 

In summary, we used neutron scattering to de- 
termine the origin of quantum paramagnetism in 
Cu(Quinoxaline)Br2. The methods of detailed balance 
correction, phonon subtraction, and multiple scattering 
correction were described and used to extract the rela- 
tively weak magnetic signal from the strong nuclear inco- 
herent scattering for quantitative comparison to theoret- 
ical models. The observed spin gap A « 1.9 meV is con- 
sistent with magnetic susceptibility and high-field mag- 
netization measurement.^^ The neutron scattering data 
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FIG. 10: (a) One-triplon contribution to the powder-averaged 
dynamic spin correlation function for a spin ladder as cal- 
culated by the Continuous Unitary Transformation. Lattice 
spacings are those appropriate for Cu(quinoxaline)Br2 and 
the exchange constants were derived by fitting to neutron 
scattering data as described in the text, (b) The comparison 
of triplon dispersion relation for the best fit between CUT 
(solid line) and SMA (dashed line) ladder models. 



are well described by the Continuous Unitary Transfor- 
mation theory of a spin- 1/2 ladder along the b-direction 
of Cu(Quinoxaline)Br2 with exchange constants J\\ — 
1.98 meV and J± = 3.3 meV. There is also evidence 
of inter-ladder dispersion with a bandwidth of approx- 
imately 0.2 meV. A c-axis alternating chain model was 
shown to be inconsistent with the data. While there may 
exist other models that are consistent with the present 
powder data, our results strongly support the interpreta- 
tion of Cu(Quinoxaline)Br2 as a quantum spin-1/2 lad- 
der with significant inter-rung coupling (J||/Jj_ = 0.6). 
Efforts are now under way to test this conclusion and ex- 




Q (A) 



FIG. 11: Comparison of inelastic magnetic neutron scatter- 
ing I(Q,u)) from Cu(Quinoxaline)Br2 at T=1.4 K (c) to (a) 
the One- Triplon scattering derived from the Continuous Uni- 
tary Transformation (b) a quasi-one-dimensional spin ladder 
model in the Single-Mode Approximation. The models were 
complemented with incoherent elastic double scattering, ac- 
cording to Eq. ||HJ with T = 0.40 and an overall additive con- 
stant was fit to account for any discrepancies in background 
subtraction. 



paramagnet through inelastic neutron scattering from an 
assembly of small crystalline samples. 
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